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Abstract. Motivated by recent numerous works on the interplay among various measures of quantum 
correlations, we aim to investigate the relationship between the violation of Clauser-Horne-Shimony-Holt 
(CHSH) Bell inequality and geometric measure of quantum discord for two-qubit systems. Exact lower and 
upper bounds of Bell violation versus geometric discord are obtained for a specific and significant class of 
states, Bell diagonal states, and the respective states which suffice those bounds are also characterized. The 
dynamical behavior of these two quantifiers is carefully analyzed in the presence of decoherence, including 
Markovian, non-Markovian, and non-back-action quantum environments. The results suggest that Bell 
violation is closely related to geometric discord, like its relationship with other entanglement monotones. 
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PACS. 03.65.Ud Entanglement and quantum nonlocality 
statistical methods - 03.67.-a Quantum information 



03.65.Yz Decoherence; open systems; quantum 



1 Introduction 

Quantum nonlocality, revealed by the violation of Bell- 
type inequalities [1 j, has already been recognized as the 
fundamental resource of present-day quantum information 
science, and quite recently, Bell violations (especially the 
violation of CHSH inequality [2]) find expanding applica- 
tions in new quantum information tasks, such as device- 
independent random number generation (DIRNG) [3] and 
quantum key distribution (DIQKD) [4,5,6 . The general 
relations between the violation of CHSH inequality and 
entropy, purity (mixedness), or several entanglement wit- 
nesses have been widely investigated [7,8,9 . So far, our 
knowledge is far from being the simple fact that violating 
a Bell inequality implies some amount of entanglement. In 
particular, for pure states, the presence of entanglement 
guarantees violation of a Bell inequality (Gisin's Theo- 
rem) [10 j. However, for mixed stares the situation becomes 
more complicated, and a great deal of endeavor has been 
devoted to this subject [11, 12, 13[ll4j. 

On the other hand, quantum discord, proposed as a 
measure of quantumness of correlations by Ollivier and 
Zurek [15 , and independently by Henderson and Vedral 
[16], has also received much attention from both theoret- 
ical [17] and experimental [18] aspects. Quantum discord 
is aiming at capturing the nonclassical part of correla- 
tions, which includes entanglement, playing a crucial role 
in some quantum information processing, especially in the 
absence of entanglement [19] • Nevertheless, the analytical 
expression for quantum discord is only available for Bell- 
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diagonal states [20] and a certain class of X-structured 
states [21 . The most recent work [22] shows that a closed 
expression for quatum discord of arbitrary two-qubit states 
cannot be obtained. Dramatically this computational dif- 
ficulty in turn motivated the proposals of alterative mea- 
sures of quantum correlations, one of which is introduced 
by Dakic et al. using the Hilbert-Schmidt norm, called 
" geometric measure of discord" [23] . Employing this defi- 
nition, an analytical formula was obtained for general two- 
qubit states in contrast to the original version of discord. 

As described above, the relations among nonlocality 
(Bell violation), mixedness, entanglement have been in- 
tensively studied. Intuitively, Bell violations should also 
show some connection with quantum discord in some form 
of mathematical expressions, since both quantities are put 
forward to characterize some kind of quantum correlations 
contained in bipartite systems from different perspectives. 
However, the related work can be rarely found in the lit- 
erature [24,25 . In the present work, we are seeking to 
address the relationship between geometric quantum dis- 
cord and nolocality, quantified by the maximal violation of 
the Bell inequality in its CHSH form. From the geometric 
point of view, we have analytically proved that the viola- 
tion of the CHSH inequality for given geometric discord 
is bounded by the relation ^\[Dg < B < 2y/l + 2D G for 
general Bell diagonal states, which is in agreement with 
the numerical result presented in Ref. [23] . To gain further 
insights, we have explicitly illustrated the dynamic pic- 
ture of Bell violation and geometric discord under several 
types of decoherence and our results indicate that these 
two quantities were closely associated with each other. 
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The remainder of the paper is arranged as follows. In 
Sec. II, we give a brief review on the notation and defini- 
tions that will be exploited throughout the paper. In Sec. 
Ill, we deduce the exact lower and upper bounds of Bell 
violation for given geometric discord. In Sec. IV, we study 
the dynamical behavior of these two quantifiers under sev- 
eral different kinds of decoherence processes. Finally, Sec. 
V is devoted to the discussion and conclusion. 



2 Preliminaries 

As a starting point, we introduce the notation and defi- 
nitions that will be employed in the remainder of the pa- 
per. Consider the two-qubit system on the Hilbert space 
%ab = C 2 (g) C 2 , the Bell operator corresponding to Bell- 
CHSH inequality can be formulated in the following form 

Bchsh = a • a ® (b + b') • a + a' • a <8> (b - b') ■ tr, (1) 

where a, a', 6, b f are the unit vectors in R 3 , and a = 
02? 03) with 01, 02, 03 being the Pauli matrices. One 
can write an arbitrary two-qubit state in the Bloch de- 
composition 

1 3 
p=-(I®I + x-a®I + I®y-a + ^ UjCJi <g> 0j){2) 

where cc, y are vectors in R 3 , and Uj = Tr(pai <S> <Jj) are 
the real components of the correlation matrix T. Then the 
well-known CHSH inequality is expressed as 

B = \(Bchsh) p \ = \Tr(pB C HSH)\ < 2. (3) 

In Ref. [26] the Horodecki family presented the necessary 
and sufficient condition for violating the CHSH inequality 
by an arbitrary two-qubit state. 

Lemma 1: (Horodecki 's Theorem f2hlj ) For any density 
matrix the maximal violation of the CHSH inequality 
ma,xTr(pBcHSH) is given by 2^/m(p), and the inequality 
is violated by some choice of a, a' , b, b' if and only 
if m(p) > 1, where m(p) = maxi<j(wj + Uj) and Ui, i = 
1,2,3 are the eigenvalues of U = T T T (T T denotes the 
transposition ofT). 

Meanwhile, geometric measure of quantum discord is 
defined as [23] 

V G (p) :=min||p- X || 2 , (4) 

where Q denotes the set of zero-discord states and \\p — 
x|| 2 = Tr(p-x) 2 is the square of Hilbert-Schmidt norm of 
Hermitian operators. In the two-qubit case, the geometric 
measure of quantum discord of Eq. (|2]) can be evaluated 
as [23] 

V G (p) = \(\\xf + \\T\\-k max ). (5) 

where k max is the largest eigenvalue of matrix K = xx T + 
TT T . 



3 Bell violation versus geometric measure of 
quantum discord 

In this section, we restrict our attention to a specific class 
of states, Bell diagonal states, not only because of analyt- 
ical simplicity but also their important role in entangle- 
ment concentration [27] , local filtering operations [7] , and 
security proof of DIQKD [6] et al. In Bloch representation, 
Bell diagonal states can be expressed as 

1 3 

pBell =-(! ® I + y^CiQ-j (g)0i) (6) 

=A+|^+)(^ + | + A+|^+)(^+| 

+ Ai|#->(#-| + ^|!P-><!?-| (7) 

where |^±) = _^ ( | l) ± |10», |$±> = ^(|00) db |11)) 
are the four Bell states, and q are real constants fulfilling 
certain constraints such that pBell is positive semi- definite, 
that is 



0< = 




- Ci 


- C 2 


- cs) < 1, 


(8a) 


0< = 




- Ci 


+ C 2 


+ C3) < 1, 


(8b) 


0<A+ = 




+ Ci 


- c 2 


+ C3) < 1, 


(8c) 


0<A+ = 




+ Ci 


+ C 2 


- cs) < 1, 


(8d) 



Utilizing these conditions, we present our Lemma 2 as a 
geometric representation of Bell diagonal states. 

Lemma 2: (Geometrical representation) If a Bell diag- 
onal state is physical, then the following restrictions must 
be satisfied: |q| < 1 for i = 1, 2, 3 and c 2 — c 2 — c| > —1 
with i 7^ j 7^ k. 

Proof: Form inequalities (|8aj) + (JHEJ) and (JHcj) + (jH2j), 
we can easily find 

-3<ci<l, -l<ci<3 =^|ci|<l, (9) 

Besides, from (|5oj) x (|8d|) and ([85j) x (jHcj), we have 

(l-cs) 2 -^!^^) 2 >0, (10a) 
(l+C3) 2 -( Cl -c 2 ) 2 >0, (10b) 

Consequently from (|lQa|) + JlOE]), we can obtain 

4-ci-4>-i. (11) 

Similarly, other constraints can be verified: | C2 1 < 1, | C3 1 < 
1 and c 2 — c 2 — c 2 > — 1 , c 2 — c\ — c 2 > — 1 due to symmetry 
property of the 3-tuple (ci, 02,03). 

In fact, in Ref [28] the Horodecki family already ad- 
vanced the geometrical representation of Bell diagonal states, 
which depicts a tetrahedron with vertices (1,1,— 1),(— 1,-1, 
(1,-1,1) and (—1,1,1). Recently, this geometric interpre- 
tation has been expanded in some other context [29] . Be- 
fore presenting our main result, we first consider the re- 
lationship between B and Dq for pure states since they 
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possess the maximum nonlocality in the sense that "entan- 
glement" just indicates "Bell violations" in sharp contrast 
to mixed states [10]. Up to local unitary equivalence, any 
pure two-qubit state can be written as 



|<9) = cos 6 |00) +sin<9 |11), 



(12) 



To calculate Bell violation and geometric discord of pure 
states, it is convenient to decompose p = \0)(6\ in Bloch 
form 



p =i(7 + cos 2(9 a z <8> I + cos 26 I <8> a z 



with the correlation components 

c\ — sin 2(9, C2 = — sin 2$, C3 = 1. 
It is straightforward to obtain 



B =2\/l + sin 2 2(9, 

D G = ^(cl + c 2 2 + c 2 3 + cos 2 20 
,2 2 2 1 _ 2 , 



1 



• max{c l5 c 2 , c 3 + cos 20}) 
sin 2 2(9, 



(15) 



(16) 



2)The upper bound. Without loss of any generality, 
one can assume cf > c| > c§, and then check the quantity 

A = 1 + 2£> G - ^B 2 



= l+2( C 2 + C 3)-(Cl+^) 

= l + i(ci-2c 2 -ci) 

> 1 + -(— 1 — 1) (using Lemma 2) 
= 0. 



(20) 



Note that in the proof we have employed the constraints 
(13) ~ c i ^ — 1 an d c i — c i — c l > _ 1- Inequality ([20]) is just 
equivalent to B < 2y/l + 2Dg, which is just the upper 
bound as we expect. 



(14) 




where c = max(ci, C2, C3). Obviously, we arrive at the re- 
lation 



B = 2y/l + 2D G . 



(17) 



Fig. 1. (Color online) The region of possible maximal Bell vio- 
lation versus geometric discord for Bell diagonal states. The or- 
ange solid curve represents Werner states while purple dashed 
curve stands for rank 1 or 2 Bell diagonal states. 



As described above, it is reasonable to conjecture that this 
relationship holds as upper bound for all Bell diagonal 
states. Now we can formulate the main theorem of this 
paper. 

Theorem: For Bell diagonal states, the violation of the 
CHSH inequality for given geometric discord is bounded 
by the relation: < B < 2^Jl + 2D G (see Fig. Qp. 

Proof: l)The lower bound. Using Lemma 1, for Bell 
diagonal states ©, CHSH violation is simply given by 

B = 2y/m(p), 
m(p) = max{c 2 + c 2 ,, c\ + c 2 , c\ + c 2 }, (18) 

On the other hand, geometric measure of discord can be 
evaluated as [23] 

D G = ^(c 2 + c 2 + c 2 -c 2 ) 

= i min{c 2 + cl, c\ + c\, c\ + c 2 }, (19) 

It is easy to see that B > A^/D G . 



Furthermore, we notice that the lower bound is satis- 
fied by Werner states [11] 

Pw = c\9-)(9-\ + {\-c)^ ce [0,1] (21) 

Meanwhile, the upper bound is attained for rank 1 or 2 
Bell diagonal states (c\ = ±1 and C2 = TC3) 

p = i±£S|#±)<*±| + I^£i |*±><*±|, (22) 

It is interesting to note that these states ([22]) are just 
states which play an essential part in sudden transition 
between classical and quantum decoherence [30]. In Ref. 
[7], Verstraete and Wolf presented the maximal violation 
and minimal violation of the CHSH ine quality for given 
concurrence C: max(2, 2^2(7) < B < 2\/\ + C 2 , and spe- 
cially for Bell diagonal states, the region of possible vio- 
lations becomes 

2 ^ 2 f +1) < B < 2VT^. (23) 

where the lower and upper bounds are sharp for Werner 
states and rank 2 Bell diagonal states respectively, which 
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is remarkably consistent with our results. Note that for 
Werner states ([2T]h the concurrence C(p) = max{0, (3c — 
l)/2}, and correspondingly geometric discord Dq = (c 2 + 
c 2 + c 2 — c 2 )/4 = c 2 /2, such that we have the relation 
4y^ = 2V2(2C(p) + l)/3 (if c > 1/3). 

4 Dynamical behavior 

In this section, we investigate the connection between Bell 
violation and geometric measure of discord in dynamical 
context. Quantum discord, as a measure of correlations 
of bipartite systems, will inevitably decohere due to the 
interaction between systems and environment. Both the- 
ory and experiment have been put forward to examine 
the dynamics of quantum discord under Markovian [3Tjl32] 
and no n- Markovian [33 ,34 environments. It is remarkable 
that quantum discord always shows a certain degree of ro- 
bustness against decoherence where in same cases entan- 
glement may display the phenomenon of " sudden death" . 
Besides, the behaviors of geometric version of quantum 
discord under various decoherence have also been stud- 
ied in the literature [35 ,36 . Recently, a number of papers 
have been devoted to the analysis of effects of environmen- 
tal decoherence on the Bell-CHSH inequality violation for 
qubit systems as well [37] . However, these two quantities 
have rarely been considered under the same dynamic con- 
ditions. In this work we focus on the comparison of their 
performance under several typical decoherence processes 
so as to verify whether they are related to each other in a 
particular way. 

4.1 Markovian noise 

We first take the phase damping (or equivalently phase 
flip) channel for example. Following the Kraus operator 
approach, the evolved state of a two-qubit system pab 
under local environments can be modeled in the Kraus 
representation [38] 

s(p) = J2K iJ p(0)Kl j , (24) 

where Kij = Kf <g> Kf are Kraus operators, satisfying 

KjjKij = J if the quantum operation is trace-preserving, 
and the operators denote the single-qubit decoher- 

ence effects. For phase damping channel, the Kraus oper- 
ators are given by 

where the parametrized time p is responsible for a wide 
range of physical phenomena [39 . Here, we restrict our 
consideration to extended Werner-like (EWL) two-qubit 
initial states 

p* =r\V){V\+ l -^h, (26) 



where r quantifies the purity of EWL states, and when 
r = 1 EWL states reduce to Bell-like pure states 

|<P> = a|01> + /?|10), |#> = a|00> + (27) 

with real parameters a and j3 defining the degree of en- 
tanglement of pure part of EWL states and a 2 + f3 2 = 1. 
Employing concurrence [40] as the measure of entangle- 
ment, the initial entanglement for both EWL states of 
Eq. ([26]) is the same, given by 

c & = c * = 2 max { , (a/3 + l/4)r - 1/4}, (28) 

Therefore it is easy to find that the initial EWL states 
are entangled only if r > (1 + 4a/3) _1 > 1/3. We plot the 
initial Bell violation and geometric discord in Fig. [2j 




Fig. 2. (Color online) Bell violation B and geometric discord 
Dg for initial EWL states as a function of parameters a 2 and 
r. 




Fig. 3. (Color online) Bell violation B and geometric discord 
Dq for decohered EWL states as a function of various param- 
eters : (a)(b) a 2 = (3 2 = 1/2 (Werner-like states); (c)(d) r = 1 
(Bell- like pure states). 

In Appendix (|6.2|) . we obtain the explicit expressions 
of B and Dq for EWL states under phase damping chan- 
nel. In Fig.[3j we show the effect of phase damping on EWL 
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states for different parameters. For Werner-like states (a 2 = 
ft 2 = 1/2), the "violation sudden death" |5HH2] clearly 
arises in contrast to the well-known phenomenon "entan- 
glement sudden death" which makes its appearance in a 
wide range of physical situations [43 j. On the contrary, 
geometric measure of discord shows its robustness against 
decoherence behaving similarly to the original version of 
discord [32]. For Bell-like states (r = 1), we note that 
the dynamics of B and Dq are almost the same, and the 
maintenance of Bell nonlocality (B > 2) is more sensitive 
to the initial purity r than a 2 and symmetric with respect 
to a 2 = 1/2. Moreover, it is remarkable to see that the re- 
lation B = 2\Jr 2 + 2Dq always holds for both initial and 
decohered EWL states (see appendix for more details). 



4.2 non-Markovian noise 

In this subsection, we analyze the non-Markovian dynam- 
ics of Bell violation and geometric discord concerning two- 
qubit system, each locally and independently interacting 
with a zero-temperature resevoir. The exact solution of 
single-qubit dynamics can be found in Ref. [44], where 
the physical process is in fact modeled as the amplitude 
decay channel 



p(t) 



_ ( Pn(0)Pt 



Pw(0)VP t 



^oi(0)VPi poo(0) + M0)(i - Pt) 
and the function Pt takes the form 

2 



-At 



. dt A , . dt 
cos( _ ) + _ sm (_) 



(29) 



(30) 



where d = V2rX — A 2 , and the parameter A denotes the 
spectral width of the coupling, r represents the decay 
rate of the excited state in the Markovian limit of a fiat 
spectrum. 

Exploiting the calculations in Appendix (j6.3ft , it is con- 
venient to implement any dynamical simulation with re- 
spect to parameters rt, r, or a 2 . For instance, we can 
easily retrieve and extend the results in Ref. [JT]. From 
Fig. |4j we see that "violation sudden death" and "viola- 
tion sudden birth" or so called "Bell islands" [41] occurs 
for both initial Werner-like states and Bell-like states. For 
Werner- like states, the dynamics of p® and p^ is just the 
same; however, for Bell-like states, the dynamics of 
is symmetric with respect to a 2 = 1/2 while pF is not. 
Note that the sudden death and revival of Bell violation 
is closely related to the memory effect of non-Markovian 
environment. If the non-Markovian effect is not sufficiently 
strong, for example, XjF is larger than 10 -2 , thus the re- 
vival can never appears. 

Futhermore, for comparison, we plot the time evolu- 
tion of the Bell violation B — 2 and normalized geometric 
discord 2Dq (the maximum value of Dq is 1/2 for two- 
qubit states) for the same fixed value of parameters in Fig. 
[5] From Ref. [44 , we know that the function P t only has 
discrete zeros at t n = 2[n7r — arctan(d/A)]/d where n is in- 
teger. The plot clearly shows the quantity B — 2 and 2Dq 




Fig. 4. (Color online) The dynamics of Bell violation 
((a)(b)) and B® ((c)(d)) in non-Markovian case as a function of 
the dimensionless quantity rt, r, or a 2 (X/T = 10 -3 ): (a)(c) 
a 2 = (3 2 = 1/2 (Werner- like states); (b)(d) r = 1 (Bell- like 
pure states). 



K W r» //\\ /'A' 




(a) 



(b) 




Fig. 5. (Color online) Bell violation B — 2 (orange solid line) 
and normalized geometric discord 2Dq (green dashed line) in 
non-Markovian case (X/T = 10 -4 ) as a function of the di- 
mensionless quantity Ft for p® ((a)(c)) and p^ ((b)(d)): (a)(b) 
a 2 1/3, r = 1 ; (c)(d) a 2 1/2, r 0.85. 



act similarly along with the function P t : simultaneously 
and periodically reach zero according to the zero points 
of P t and achieve extreme values at the maximum points 
of Pt, no matter which states we consider, p® or p^ (one 
can easily obtain the exact expressions of B — 2 and 2Dq 
from Appendix (j6.3p so as to justify all the details). These 
figures show that times of revival of Bell violation sensibly 
depend on the purity parameter r, and we see that times 
of revival reduce if r becomes smaller. However, geometric 
discord only vanishes at discrete points and still remains 
robust against non-Markovian noise like the situation in 
the Markovian case. 
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4.3 non-back-action environment 

As has been mentioned above, the occurrence of Bell vi- 
olation revivals originates from the memory effect of the 
non-Markovian quantum reservoirs, or more precisely, due 
to the back-action of qubit system on the corresponding 
environment. This back-action mechanism leads to estab- 
lishment of quantum correlations between qubits and envi- 
ronments or between the environments themselves. Here, 
we investigate the case where a pair of independent qubits 
each locally subject to the so called "system-unaffected 
environment" (SUE) [35], which means the back-action of 
system does not exist, in order to see whether the revivals 
of Bell violation or geometric discord would happen in this 
model. 

Consider a initial two-qubit system p(0) each locally 
coupled to random external fields, which can be charac- 
terized as a bistochastic quantum channel [36]. The dy- 
namical map A can be always written as [17] 

N 

A[ P (0)] = Y.PiAipMA*, (31) 

2=1 

where pi is a probability measure with U-L-^Pi = l,Pi>0 
and Ai, i = 1,2,..., TV are unitary operators. The environ- 
ment depicted in Ref. [45] can be imagined as follows: two 
classical external fields have the same frequency and am- 
plitude but each pass through a random dephaser which 
switches between the values and tt with the same prob- 
ability, and then locally act on one of the qubit pair. For 
this case the single-qubit dynamical map can be written 
as A s [ps(0)} = \ E- =1 Uf (t)p s (0)C/f \t), where Uf(t) = 
e -iHit/h foe time evolution operator with Hi = Huja z /2-\- 
g(a+e~' l( f )i -f cr_e^) and can be represented in matrix form 
in {|1), |0)} basis 

tjS (A _( cos(gt) -e-*** sm(gi)\ , , 

Ui[t) ~\e i ^sm(gt) cos(gt) )> [6Z) 

Successively, the time behavior of the two-qubit system 
Pab can be given by 

P {t) = A[p(0)] = \Y, U*(t)Uf{t) P {0)U^uf\ (33) 

A straightforward calculation then shows that the map A 
works inside the class of Bell-diagonal states, and partic- 
ularly the effect of this map on Bell states have the form 

W±)(P±\ = [1 " f{t)\\P±){P±\ + f{t)\p' T ){(3' T \, (34) 

where /3,/3' = with /3 ^ /3' and \&±) = (|01) ± 
|10»/V% |#±) = (|00)±|11))/V2, and /(*) = sm 2 (2gt)/2. 
For a general Bell-diagonal state p = ^2 a s \p(Q)\/3 3 )(/3 8 \ 
(/? = \P, <P; s = =b), the time-dependent coefficients of p(t) 
can be obtained as 

Af (*) = A±(0)[1 - /(*)] + Aj,(0)/(t), (P + /*') (35) 



0.8 r 

0.7 




Fig. 6. (Color online) The time behavior of Bell violation B — 
2 (orange solid line) and normalized geometric discord 2Dq 
(green dashed line) as a function of the dimensionless quantity 
gt for an initial Bell-diagonal states A+(0) = 0.9, A 1 "(0) =0.1. 

Though it is easy to give the explicit formula for Bell 
violation or geometric discord based on all the unspecified 
parameters A^(0), we now focus on the initial conditions 

A+(0) = 0.9, A^(0) = 0.1, and A+(0) = A^(0) = 0. Utiliz- 
ing Eq. (j35j), we have 

A+ = 0.9[l -/(*)], =0.1[1 -/(*)], 

A+ = 0.1/(t),Ai = 0.9/(t), (36) 

Note that for Bell-diagonal states in Bell basis, the maxi- 
mum possible violation reads as [7([T^] 

B = 2v / 2\/(A 1 - A 4 ) 2 + (A 2 - A 3 ) 2 , (37) 

with eigenvalues Ai > A2 > A3 > A4. In our case, A^ is 
largest and A J is smallest among the eigenvalues, due to 
f(t) G [0, 1/2]. Therefore the Bell violation is given by 

B = 2v / 2V[0-9-/W] 2 + [0.1-/(t)] 2 , (38) 
The corresponding Bloch representation is 

ci = 0.8 - 1.6/(t), c 2 = 0.8, c 3 = 2f(t) - 1, (39) 
Then the geometric discord can be cast as 

r (0.8-1.6/) 2 +0.8 2 </<01 
DG = \ (0-8-l-6/)V(l-2/)\ 7 " Q 5 ( 4 °) 

Fig. [6] displays that Bell violation B — 2 periodically 
vanishes and revives, with its revival amplitude unchanged. 
This behavior is in sharp contrast with the situation in 
the non-Markovian case, where the revival amplitude of 
Bell violation is damping gradually until the complete 
death. The maximum values of B and Dq are attained 
at t n = nir/2g (n = 1,2,...), while Dq only vanishes at 
t n = (2n — l)7r/(4g) like its original version D showed in 
Ref. [35] . It is worth noting that the evolution of the rele- 
vant quantifiers is of great dependence on the initial state 
in this environment. For instance, if all the initial eigen- 
values A^(0) (f3 = <F, s = ±) are less than 1/2, then the 
revival of Bell violation will never occur. 
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5 Conclusion and Open questions 

First, We discuss the relationship between the violation of 
CHSH-Bell inequality and geometric measure of quantum 
discord for Bell diagonal states. Exact lower and upper 
bounds of Bell violation versus geometric discord are de- 
rived and the respective states which suffice those bounds 
are also characterized: the lower and upper bounds are at- 
tained for Werner states and rank 2 Bell diagonal states 
respectively. Second, the dynamical behavior of these two 
quantifiers is analyzed in the presence of decoherence, in- 
cluding Markovian, non-Markovian, and non-back-action 
quantum environments, especially for extended Werner 
states. Under Markovian noise (phase damping channel), 
sudden death of Bell violation occurs concerning with the 
purity parameter r, but the revival never appears. More- 
over, it is interesting to find that the relation B = 2\Jr 2 + 2Dq 
always holds for both initial and decohered EWL states 
(comparing with the well-known formula B = 2^1 + 2Dq 
for pure states). In the non-Markovian case, Bell violation 
periodically vanishes and revives according to the function 
Pt with its revival amplitude damping, while geometric 
discord only vanishes at some discrete points. More im- 
portantly, we have clearly displayed that B — 2 and Dq 
simultaneously reach zero according to the zero points of 
Pt and achieve extreme values at the maximum points 
of P t . Note that the sudden death and revival of Bell 
violation is closely related to the memory effect of non- 
Markovian environment (the value of X/T). Finally, we 
consider the time evolution of B — 2 and Dq driven by 
a system-unaffected environment. Our result shows that 
Bell violation B — 2 also periodically vanishes and revives, 
however, with its revival amplitude unchanged, which is 
in sharp contrast with the situation in the non-Markovian 
case. It is worth emphasizing that the maximum values of 
B and Dq are attained at the same time, and the evolu- 
tion of these two quantifiers is of great dependence on the 
initial state in this environment. 

Yet there are still some open theoretical problems for 
further research. One is the question that whether there 
exists a strict lower or upper bound of Bell violation ver- 
sus geometric discord for more general states, we conjec- 
ture that B < 2y/l + 2Dq holds for almost all states, but 
the proof may need more technical tricks. In addition, we 
are seeking to investigate the relation between Bell vio- 
lations (non-locality) and other correlation quantities, for 
example, randomness, concerning that the Bell inequality 
violation can be used to certify the presence of genuine 
randomness [3|[48] . 
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6 Appendix 

6.1 B and Dq for X states 

In this appendix, we analytically present the maximum 
violation of Bell-CHSH inequality and geometric measure 
of discord for a more general class of two-qubit states, X- 
structured states, which are represented in the orthonor- 
mal basis {|00), |01), |10), |11)} 



Px 



( Pii 



\P41 





P22 P23 
P32 P33 




Pl4\ 




pmJ 



(41) 



where we can assume P23 — P32 and P14 = p±\. For com- 
putational simplicity, it is helpful to rewrite the state ([4T] 
in Bloch decomposition 



>/ + ra-cr 3 (g)/ + /(g)n-cr3 



3 

E< 

i=l 



(42) 



with 



ci 

C2 
C3 

rn 



= 2pi4 + 2p 2 3, 

= -2/9U + 2p 2 3, 
= Pll - P22 - P33 + P44, 
= Pll + P22 — P33 — P44, 
Pll — P22 + P33 — P44- 



(43) 



2^ 



According to the Horodecki criterion [26] , B 
with i : j = 1,2,3. The three eigenvalues U{ of U = T l T 
are 



Ui 
U2 



U3 



4(|pi 4 | + |p23|r, 
4(|P1 4 |-|P23|) 2 , 
{pll — P22 — P33 



P44) 2 



(44) 



It is easy to see that u\ is always larger than ^2, and thus 
the maximum violation for X states is 
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B = max(Bi, B 2 ), B 1 = 2y/m + u 2 , B 1 = 2y/uT 



f^3, 

(45) 



Meanwhile, geometric measure of discord can be evaluated 
as 



d g = -Ac 



4 



m 



max{cf , c|, C3 + m 2 }). 



(46) 
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6.2 Markovian noise 

In the standard product basis {|11), |10), |01), |00)}, we 
can rewrite EWL states in matrix form 



p*(0) = 



/ 


1-r 
4 








° ^ 









aj3r 










af3r 







V 











1-r 
4 / 



(47) 



/i^ + /? 2 r 



p*(0) 









y aj3r 



^ o 

1=E 



af3r \ 





(48) 



±=z + a 2 r / 



Note that EWL states are a subset of X states, and thus 
they can also be represented in the form of Eq. (|^2j) . For 
p^(0), these parameters are given as follows 



ci — 2af3r, C2 = 2af3r, cs 



-r, m = —n = (f3 2 — a 2 )r. 

(49) 



While for p^(0), we have 

c\ = 2af3r, C2 = — 2a/3r, C3 = r, m = n = (/3 2 — a 2 )r. 

(50) 

Keeping a 2 + j3 2 = 1 in mind, it is straightforward to 
obtain the same B (and Dq) for both p^(0) and p*(0) 



5 = 2 v / r 2 +4a 2 ^ 2 r 2 , L> G = 2a 2 p 2 r 2 . (51) 
As is the case for phase damping, EWL states evolve as 



/(p) 



p*(p) 



/i^ \ 

^+j3 2 r (l-p)aj3r 

(l-p)a/?r ^+a 2 r 

v o o o 



f^T+^r (l-p)o;/3r\ 
0^0 
if£ o 

\ (1 - p)a(3r ±^ + a 2 r ) 



(52) 



(53) 



Through Eq. (|24|) and ([25]) , we can find the only change 
is that c\ and C2 become = (1 — p)c\ and c' 2 = (1 —p)c2 
respectively for both p^(p) and p^(p), and other param- 
eters C3,m,n remain unchanged. Hence, 5 (and Dq) are 
still the same with respect to p®{p) and p*(p) (note that 
<p < 1) 



5 = 2Vr 2 +4(1 -p) 2 a 2 /3 2 r 2 , L> G = 2(1 - p) V/3 2 r 2 . 

(54) 



6.3 non- Markovian noise 

Following the non-Markovian model described in Ref. [33] , 
we can easily obtain the exact evolution of EWL states. 



For p<£, the decohered nonvanishing density matrix ele- 
ments are 

pf 2 (i) = + /? 2 rm + - p t ), 

pUt) = + o?r)P t + ^P t (l - Pt), 

p* 4 (t) = l-P t + 1 -^P 2 , 



(55) 



The corresponding Bloch representation is also of the form 
Eq. (021) 

c 1 (t) = c 2 (t) = 2aprP t , 

c 3 (t) = l-2P t + (l-r)P 2 , 

m(t) = (/? 2 - a 2 )rP t - (1 - P t ), 

n(t) = (a 2 - /3 2 )rP t - (1 - P t ), (56) 

Meanwhile, for p^, we obtain 

PU) = + P 2 r)Pl 

Pt2(t) = pfs(t) = l -^P + (^ + /? 2 r)P t (l - P t ), 

PUt) = 1 + + fr)P 2 t - (1 - r + 2/3 2 r)P t , 

= P*i{t) = <*PrP t , (57) 
The corresponding Bloch decomposition is given by 
ci(t) = -c 2 (t) = 2a/3rP t , 

c 3 (t) = 1 + (^ + p 2 r)4P t (l - P t ) - (1 - r)P t , 



m(t) = n(i) 



?2 „2 



a>P t - (1 - P t ), 



(58) 



Exploiting the above formulas, we can conveniently carry 
out any dynamical simulation with respect to all these 
parameters. 
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